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Abstract 

In this paper, we study the control of the linear heat equation with a space and time 
dependent coefficient function by the Dirichlet and Neumann boundary control laws. 
This equation models the heat diffusion and space, time dependent heat generation in a 
one dimensional rod. Without control, the system is unstable if the coefficient function 
is positive and large. With boundary control based on the state feedback, we show that 
for the time analytic coefficient a(x,t), the exponential stability of the system at any 
rate can be achieved. It is further shown that both the control of the Dirichlet and 
Neumann boundary value systems can be stabilized using this method. In doing this, 
the control kernels are explicitly calculated as series of approximation and they are used 
in the simulations. The numerical simulation confirmed the theoretical arguments and 
the controllability of the system. The possible future works are discussed at the end. 
This is the continuation of the recent work of Liu [SIAM J. Cont. Optim. vol. 42, pp. 
1033-1043] and Smyshlyaev and Krstic [accepted by Automatica]. 
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1 Introduction 



In this paper, we continue the study of feedback control of an unstable heat equation 



This equation arises from the heat generation and diffusion of a rod in which the heat gen- 
eration is only only sensitive to the position x but also the time t. A simple version of this 
feedback control problem where a = fi, a constant, was first addressed by Boskovic, Krstic 
and Liu in pQ and the instability when fi is large was shown there. Later, Balogh and Krstic 
[21 E! studied the case where a = a(x), a function dependent on x. Since they used a backstep- 
ping approximation method, the same problem was reconsidered and solved by Liu [IJ using 
a transformation skill. More recently, the generalization to a = a(t) was explicitly considered 
by Smyshlyaev and Krstic j^J. In our paper, we generalize it to the case of a = a(x,t). 

For the theories that use boundary feedback control to stabilize the parabolic equations, 
one can refer to j^j and |3] and references there. In comparison of the existing literature, the 
novelty of this paper is the "explicit" construction of the feedback laws of a = a(x,t) case 
for some a(x, t) and the proof that both the kernel for the Dirichlet and Neumann boundary 
control problem exist. 

The paper is organized as follows. Section El is devoted to the stabilization of Dirichlet 
boundary value problems of a = a(x, t) case, Section |3] to the stabilization of Neumann 
boundary value problems. In Section HI the simulation of a(x, t) = x{bt + c) is studied and 
the validity of the control are shown. In section 5, we give some remarks of interesting 
phenomena, which may be helpful for the future works. 

In what follows, we denote the region (0, l) x by X, (0, l)t by T, < y < x < 1 by Q and 
X x T by XT and x T by QT and so on. We also denote by H s : XT the usual Sobolev 
space (see, e.g., |H1IID1) for any s G R. For s > 0, Hq-.XT denotes the completion of :XT 
in H S :XT, where C^°:XT denotes the space of all infinitely differentiable functions on XT 
with compact support in XT. We denote by || • || the norm of I? :XT with respect to space 
variable x. C n : XT denotes the space of all n times continuously differentiable functions on 



u t (x, t) = u xx (x, t) + au(x, t) in (0, 1). 



(1.1) 



XT. 



2 Dirichlet Boundary Conditions 



It is well known that the Dirichlet boundary value problem 




(2.1) 



is unstable if function a is positive and large 
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To design a boundary feedback law to stabilize the above system for function a(x,t), we 
consider the problem 

k xx (x, y, t) - k yv (x, y, t) - k t (x, y, t) = (a(y, t) + X)k(x, y, t), in QT 
k(x,0,t) = 0, in XT (2.2) 

k x (x, x, t) + k y (x, x, t) + 4z{k{x, x, t)) = a(x, t) + A, in XT 

where A is any positive constant. Why we consider this problem will becomes clear in 
Lemma \'2. 21 For the moment, let us assume this problem has a unique solution k for some a. 
Using the solution k, we then obtain a Dirichlet boundary feedback law 

t) = - I fc(l, y, t)u(y, t)dy in T (2.3) 
Jo 

and Neumann boundary feedback law 

u x {l, t) = -k{l, t) - [ k x (l, y, t)u(y, t)dy in T. (2.4) 

With one of the boundary feedback laws, the system 

u t (x, t) = u xx (x, t) + a(x, t)u(x, t) in XT, 

w(0,t) = inT, (2.5) 

u(x, 0) = u°(x) in T 

is exponentially stable. In this controlled system, the left hand end of a rod is insulated while 
the temperature or the heat flux at the other end is adjusted according to the measurement 
of fc-weighted averaged temperature over the whole rod. Physically, if the destabilizing heat 
is generating inside the rod, then we cool the right end of the rod so that it is not overheated. 
To state this result, we also introduce the compatible conditions for the initial data 

u°(0) = 0, u°(l,t) = -[ k(l,y,t)u°(y)dy (2.6) 



«°(0) = 0, u° x (l,t) = -k{l,l)u°(l,t) - k x (l,y,t)u°(y)dy (2.7) 



Theorem 2.1. Assume that A > is any positive constant and a(x,t) is continuous in x 
and analytic in t uniformly in x. For arbitrary initial data u°(x) G H 1 X with compatible 
condition (|2.6J) or (12. 7j) . equation ()2.5|) with either (|2.3J) or (12. 4j) has a unique solution that 
satisfies 

\\u{t)\\ H i < M\\u°\\ H ie- xt , t G (0, to), to < 1, (2.8) 
where M is a positive constant independent of u°. 

The idea of proving the theorem is to carefully construct a transformation 

/•X 

w(x,t) = u(x,t) + / k(x,y,t)u(y,t)dy 
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to convert the system (J2.5)) with either (J2.3)) or (J2.4)) into the exponentially stable system 

Wt = w xx — \w in XT, 

w(0,t) = w(l,t) = inT, (2.9) 
w(x, 0) = w°(x) in X, 

or 

w t = w xx — Xw in XT, 

w(0,t) = w x (l,t) = inT, (2.10) 

0) = u>°(x) in X, 

where w°(x) = u°(x) + J* k(x,y,t)u°(y)dy. This will be done by the following lemmas. The 
first one of these lemmas is due to Colton |H] . We conveniently rewrite its statement and 
quote the proof here for reference. 

Lemma 2.1. (Colton) Suppose that a(x,t) is continuous in x and analytic in t uniformly 
in x. Then problem ()2.2|) has a unique solution which is twice continuously differentiable in 
x, y. 

Proof. Using the substitutions 

x = t + V, V = €-ri, a(y, t) = a(£ -7},t), G(£, r], t) = k(x, y, t), 
the problem (|2.2j) is transformed to 

G^, 7i, t) = (a (f - 7], t) + A + |) V,t) 0<77<e<l,ter 

t) = 0<^<l,teT 
|(G(£,0,t) = ± (a(f,t) + A) < e < l,t 6 T. 



Integrating over £ and 77, we can find 

G{£, V,t) = \ J (a(r, + A) dr + £ jH ( 



a (r - s,t) + A + — ) G(r, s, t)rfsrfr . (2.11) 
at 



Then Colton showed that the above equation has a solution using the "method of domi- 
nants" . This method works as follows. If we are given two series 

00 00 
Si = $>i n f\ S 2 = 5> 2 X, t e (0,1), (2.12) 

n=l n=l 

where a 2n > 0, then we say S 2 dominates Si if |a ln | < a 2n , n — 1, 2, 3, • • • , and write Si <C 5*2. 
It can be easily checked that 

ifSi<S 2 , then|Si|<S 2 , (2.13) 

and ^ « Si « S 2 (l - t)" 1 ; (2.14) 
at ot 

if 5i < S 2 , S 2 < S 3 , then Si < S 3 ; (2.15) 

if Si < S 2 , S 3 < S 4 , then Si + S 2 < S 3 + S 4 . (2. 16) 
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Moreover, the property of "dominant" can also be kept if the integrals of the two series are 
not with respect to t but other variables, that is 

rb rb 

ifSi<S 2 , then / Sxdx^ J S 2 dx. (2.17) 

J a J a 

On the other hand, if a function / is analytic in t e (0, 1), then there exist a positive constant 
C such that / < C(l - t)' x . 

Using this method in our problem, it can be shown that equation (j2.11|) has a unique 
twice continuously different iable solution if a is analytic in t. In fact, since a(x,t) is analytic, 
we can let C be a positive constant such that we have 

a{£-V>t) +A<C(l-t)- 1 . 

By the properties of dominant (|2.13|) - ([2.17]h it can be shown by induction that for the fol- 
lowing series 

oo 

G(£,r,,t)=Y l Gn(Z,ri,t) (2-18) 



n=0 



where 



1 r 

G (Z,V,t) = 7,j (a(r,t) + X)dr 
G n +i(€,V,t) = f f Q (^ a ( T - s ^t) + ^ + G n (r,s,t)dsdr 



we have 



and hence 



nn/~<n tn„n 

G n «^|^(l-t)- (n+1) 



nl 



\G n \ < ) V (l-*r (w+1) • (2.19) 
nl 

Thus the series ()2.18|) converges absolutely and uniformly. Moreover, since a is C 1 , G is 
clearly C 2 in x, y. 

□ 



Remark 2.1. The proof of Lemma [2.11 provides a numeric computation scheme of successive 
approximation to compute the kernel function k in our feedback laws (J2.3)) and (J2.4|) . This 
makes the feedback laws (|2.3J) and ()2.4j) very useful in real problems. 

Lemma 2.2. Let k(x,y,t) be the solution of problem (j2.2j) and define the linear bounded 
operator K : H*: XT -> W : XT (i = 0, 1, 2) where T denotes (0, t„) (*o < 1) by 



w(x,t) = (Ku)(x,t) = u(x,t) + / k(x,y,t)u(y,t)dy, for u E : XT . (2.20) 
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1. K has a linear bounded inverse K 1 : H l :XT — > H l :XT (i = 0, 1,2). 

2. K converts the system (J2~3|) and (jZSJ) and system and (J2~5|l mto (J2~9|) and pTTT)|) . 

respectwe/?/. 

Proof. To prove 1.: that ()2.20|) has a bounded inverse, we set 

v(x,t) = / k{x,y,t)u(y,t)dy 
Jo 

and then 

w(x, t) = u(x, t) + v(x, t). 

Hence we have 



v(x,t) = I k(x,y,t)[w(y,t) -v(y,t)]dy 

px 

k(x,y,t)w(y,t)dy - / k(x,y,t)v{y,t)dy. (2.21) 
Jo 

To show that this equation has a unique continuous solution, we set 

PX 

v (x,t) = / k(x,y,t)w(y,t)dy, 
Jo 

px 

v n {x,t) = - k(x,y,t)v n ^ 1 (y,t)dy 



Though from 1)2.19)1 we see that k(x,y,t) is divergent with t approaches 1~, it is still clear 
that under the requirement t e (0,t ), the k(x,y,t) is still bounded in QT . Thus we can 
denote M = sup \k(x,y,t)\. Then one can see 

0<y<x<l,t£T 

\v {x,t)\ < M\\w\\, 

\vi{x,t)\ < M 2 \\w\\x, 

. .. M 3 \\w\\ 2 
\v 2 (x,t)\ < — - } — x , 

and by induction, 

v n (x,t) < ^x n . 

nl 

These estimates show that the series 

v(x,t) = ^2v n (x,t) 

n=0 

converges absolutely and uniformly in XT and that its sum is a continuous solution of 
equation ()2.21j) . Moreover, there exists a constant C > such that 

||u||<C|H|. (2.22) 
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This implies that there exists a bounded linear operator $ : L 2 : XT — > L 2 : XT such that 

v(x, t) = (&w)(x,t) 

and then 

u(x,t) =w{x,t) -v(x,t) = {{I - $)w)(x,t) = {K- 1 w){x,t). (2.23) 
It is clear that K' 1 : L 2 : XT L 2 : XT is bounded. To show that K' 1 : H 1 : XT H 1 : 



XTq is bounded, we take derivative in (j2.21|) and obtain 

v x (x,t) 

PX PX 

= k(x,x,t)w(x,t) + / k x (x,y,t)w(y,t)dy - k(x,x,t)v(x,t) - / k x (x,y,t)v(y,t)dy, 
Jo Jo 

which, combining ()2.22|) . implies that there exists constant C > such that 

\\v x \\ < C||w|| 

and then by (l2~2l?|) 

INI/p < II^IIh 1 + IMIh 1 < cII^IIh 1 - 

By analogy, we can show that K^ 1 : H 2 :XTq — > H 2 :XTq is bounded. 

To prove 2.: that the transformation ()2.2()j) converts the system ()2.3j) and (|2.5jl and system 
(12. 4j) and (|2.5|l into ()2.9j) and (|2.1()jl . respectively, we compute as follows. 

PX PX 

w t {x,t) = u t {x,t) + k(x,y,t)ut(y,t)dy+ k t (x,y,t)u(y,t)dy 
Jo Jo 

PX PX 

= u t (x,t)+ / k(x,y,t)[uyy(y,t) +a(y,t)u(y,t))dy+ / k t (x,y,t)u(y,t)dy 
Jo Jo 

= Ut(x, t) + k(x, x, t)u x (x, t) — k(x, 0, t)u x (0, t) 
— k y (x, x, t)u(x, t) + k y (x, 0, t)u(0, t) 

PX 

+ / [kyy(x,y,t) + k t (x,y,t) + k(x,y,t)a(y,t)]u(y,t)dy, (2.24) 
Jo 



w x {x 1 1) = u x {x, t) + k(x, x, t)u(x, t) + / k x (x,y,t)u(y,t)dy, (2.25) 



o 



w xx (x, t) = u xx (x, t) + -^(k( x i X ' t)) U ( X i *) + H X i X ' t)) u x(x, 1) 



X 



+ k x (x,x,t)u(x,t) + k xx (x,y,t)u(y,t)dy. (2.26) 



o 
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It then follows from (Q and ((231) that 
w t - w xx + Xw 
= u t (x, t) + k(x, x, t)u x (x, t) — k(x, 0, t)u x (0, t) 
— k y (x, x, t)u(x, t) + k y (x, 0, t)u(0, t) 

PX 

+ / [k yy (x,y,t) + kt(x,y,t) + k(x,y,t)a(y,t)]u(y,t)dy 



u xx (x,t) - —(k(x,x,t))u(x,t) - k(x,x,t))u x (x,t) 
ax 



- k x (x, x, t)u(x, t) - k xx (x, y, t)u(y, t)dy 
Jo 

+ Xu(x, t) + X / k(x,y,t)u(y,t)dy 



Jo 

a(x, t) — k x (x, x, t) — k y (x, x, t) — —(k(x, x, t)) + X^j u(x, t) 
+ k y (x, 0, t)u(0, t) - k(x, 0, t)u x (0, i) 

PX 

+ / [kyy(x, y, t) + k t (x, y, t) - k xx (x, y, t) + (a(y, t) + X)k(x, y, t)] u(y, t)dy 
Jo 

= 0. (2.27) 

By the boundary condition of (|2.5|) . we deduce that w(0,t) = 0. Using feedback law (|2.3|) or 
(l2~l) . we obtain 

w{l,t) =u(l,t)+ [ k{l,y,t)u(y,t)dy=0, 



or 



w x (l,t) =u x (l,t) + k(l,l)u(l,t)+ / k x (l,y,t)u(y,t)dy = 0. 

Jo 

□ 

Remark 2.2. We require t G (0, to)(^o < 1) i n this lemma since k is not bounded when 
t approaches 1~. However, this requirement is tolerable since in practice, we only require 
the system to be stable in a time interval t 6 (0, T ) and this can be rescaled to (0, t ) by 
t — > t/T(T > To) the heat equation. Moreover, T can be choose as big as we want. 

Now we are ready to prove Theorem 12.11 

Proof of Theorem \2. 11 We first note that problem ()2.5j) with either (|2.3J) or ()2.4j) is well 
posed since, by Lemma I2.2| they can be transformed to the problem (12. 9 J) or ()2.10j) via the 
isomorphism defined by (|2.2()jl and the problem ()2.9j) or (|2.1()|) is well posed (see, e.g., [HI 
Chap. IV]). Moreover, there exists a positive constant C > such that 

\\u(t)\\ m < C\\w(t)\\ Hh 
\\w°\\ m < C\\u°\\ m . 
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Therefore, it is sufficient to prove (|2.8|) for the solution w of (|2.9|) or (|2.1(Jj) . We do so only 
for problem ()2.9|) since the situation for problem ()2.10|) is similar. 
We define the energy 

If 1 

E(t) = - w(x,t) 2 dx. 
2 Jo 

Multiplying the first equation of ()2.9|) by w and integrating from to 1 by parts we get 



i 



1 



E(t) = w x w — / w x (x,t) dx — X w(x } t) dx 



o JO JO 

i pi 

2j„ \ / „.J„ +\2, 



which implies 



Set 



w x (x, t) dx — X w(x,t) dx 
'o Jo 
< -2XE(t), 



E(t) < E{0)e~ 2Xt , for t > 0. 



V(t) = / w x (x,tydx. 
Jo 

Multiplying the first equation of (|2.9j) by w xx and integrating from to 1 by parts we obtain 

V(t) = — 2 / w 2 xx dx + 2 A / ww xx dx 
Jo Jo 

= —2 w 2 xx dx — 2A / w 2 x dx 

Jo Jo 
< —2XV(t), 

which implies that 

V(t) < V(0)e- 2Xt . 

This shows that $TjJ$ holds. □ 

3 Neumann Boundary Conditions 

To stabilize the Neumann boundary value problem 



u t (x, t) = u xx (x, t) + a(x, t)u(x, t) in XT, 
u x (0,t) = u x (l,t) = in T, 



(3.1) 



we consider the problem 



k xx (x, y, t) - k yy (x, y, t) = (a(y, t) + X)k(x, y, t), in QT 

k y (x,0) = 0, inX (3.2) 

k x (x, x, t) + k y (x, x, t) + -^(k(x, x, £)) = a(x, t) + A, in XT 
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where A is any constant. Using the solution k, we then obtain a Dirichlet boundary feedback 
law 

u(l,t) = - / k{l,y,t)u{y,t)dy in T (3.3) 
./o 

and Neumann boundary feedback law 

u a .(l,t) = -A:(l,lHl,t)- / k x (l,y,t)u{y,t)dy in T. (3.4) 

>/ o 

With one of the boundary feedback laws, the system 

u t {x, t) = u xx (x, t) + a(x, t)u(x, t) in XT, 

u x {Q,t) = inT, (3.5) 

u(x, 0) = u°(x) in X 

is exponentially stable. To state this result, we introduce the compatible conditions for the 
initial data 

u° x (0) = 0, u°(l,t) = - [ k(l,y,t)u°(y)dy (3.6) 



o 



u° x (0) = 0, u x (l,t) = -k(l,l)u°(l,t)- k x (l,y,t)u°(y)dy (3.7) 



Theorem 3.1. Assume that A > is any positive constant and a(x, t) is continuous in x and 
analytic in t uniformly in x. For arbitrary initial data u°(x) G -£^(0, 1) with the compatible 
condition (j3.6J) or f!3.7|) . equation f!3.5|) with either (j3.3J) or (13. 4j) has a unique solution that 
satisfies 

\\u(t)\\ H i < M||n°|| H ie- A ', t e (0,t )(to < 1) 
where M is a positive constant independent ofu°. 

Proof. The proof is the same as that of Theorem l2.ll Only thing we need to do is to show that 
problem (|3.2J) has a unique solution. This is also considered by Colton [Zj and his theorem 
is quoted as Lemma 13. II below. □ 

Lemma 3.1. (Colton) Suppose that a(x,t) is continuous in x and analytic in t uniformly 
in x. Then problem ()3.2j) has a unique solution which is twice continuously differentiable in 
x, y. 

Similar to Lemma 12.21 we have 

Lemma 3.2. Let k(x,y,t) be the solution of problem (j3.2j) and define the linear bounded 
operator K : H': XT -> H* : XT (i = 0,1,2) by 



w(x, t) — (Ku)(x, t) — u(x, t) + I k(x,y,t)u(y,t)dy. 

Then, 
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1. K has a linear bounded inverse K 1 : H l :XT — > H l :XT (i = 0, 1,2). 

2. K converts the system (|3.3j) and (|3.5|) and system (|3.4|) and (|3.5|) into 



w t = w xx 


— Xw 


in XT, 


w x (0,t) = 


-- w(l,t) = 


in T, 


w(x, 0) = 


w°(x) 


in X, 


w t = w xx 


— Xw 


in XT, 


w x {0,t) = 


w x (l,t) = 


in T, 


w(x, 0) = 


w°(x) 


in X, 



respectively, where w°(x) = u°(x) + f Q k(x,y,t)u°(y)dy. 



4 Simulations 

To see how the boundary control stabilize the heat equation, we studied the case of Dirichlet 
boundary problem (j2.1)l in which 

a(x,t) = x(bt + c), (4.1) 

where b > and c > are positive constants. Both the controlled system by the Dirichlet 
boundary control ()2.3|) and uncontrolled systems are studied. In all the simulation, we use 
the initial data 

u°(x) = 10 ^ - (^x - ^ j sin(47rx) + 5 

which is arbitrarily chosen. If not stated explicitly, the following values are used in the 
simulation: A = 10, b = 200, c = 5. The x grid is set to be 100 and the time step size is set 
to be 10~ 5 but when the figures are drawn, we used larger step sizes in x and t for clarity. 

The corresponding model of the simulation can be viewed as the heat diffusion along a 
one dimensional rod, during which heat is also generated unevenly and increasingly. The 
temperature of the left end of this rod is fixed zero and another is allowed to adjust to let 
the temperature of the rod converge to zero. 

To reveal the various features and relations of the uncontrolled system, controlled system 
and the control inputs, we do the following simulations. We studied the effects of different 
6's in a(x,t) with b = 200 as shown in figure (|5.11 a) and b = 150 in figure (|5.1l c) to the 
uncontrolled system ()2.1|) . Clearly, (1) the solution to the uncontrolled system ()2.1|) is really 
divergent with time going and (2) the larger the coefficient function a(x, t) is (here the bigger 
the b is), the faster the solution diverges. The solution to the controlled system with b = 200 
is shown in Figure (|5.11 b). It is seen that the feedback control can really stabilize the system. 
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We also studied the effect of different A's in the kernel k, with A = 10 shown in figure (|5.1l b) 
and A = 40 in figure ([5. li d). Comparing the x — 1 edge of these two figures, we can find 
that the larger the A is, the faster the solution converges to zero, as expected from equation 
(|2.8jl . In all the simulations of (|5.11 a-d). the kernels are taken to be the first 3 terms in it 
series approximation ()2.18j) . that is Go + Gi + Gi- In (15. li e), we showed the difference of 
solutions to the system controlled by a kernel of first three terms ( (|5.1l b)) and the system 
controlled by a kernel of first four terms. If denoting these two solutions by m 3 (x, t) and 
Ui(x, t), in figure (|5.1l e) we showed (u 3 (x, t) —u^Xjt)), from which we can see that the high 
order approximated kernel can stabilize the system in a faster way. However, this is not a 
general conclusion for all a(x,t) since in our case, all terms of G n are positive and thus the 
more terms you use, the better approximation and better control it gives. 

5 Remarks 

As a future work, one can try to stabilize the nonlinear problem 

Ut(x, t) = u xx (x, t) + a(x, t, u(x, t))u(x, t) 

where a depends not only on x, t but also on u itself. However, we don't know how to 
prove the unique existence of the corresponding kernel. Once the unique existence of the 
corresponding kernel was shown, all the results in section (j2J) and (J3J) may be generalize to 
this case. 

Another interesting result arises from the simulation. We tried to calculate the solution 
to the controlled system (by 3 terms kernel) after t — 1, which is out of the time domain in 
our theorems. From figure (|5.1l f). one can see the system is still stable. This implies that 
the condition that a(x, t) should be analytic in t in theorems (|2.1jl and ()3.1|) is only sufficient 
but not necessary. 
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